Phase transition from a d x 2_ y 2 to d x 2_ y 2 + d xy superconductor 
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We study the phase transition from a d x 2_ y 2 to d x 2_ y 2 + d xy superconductor using the tight-binding 
model of two-dimensional cuprates. As the temperature is lowered past the critical temperature 
T c , first a d x 2_ y 2 superconducting phase is created. With further reduction of temperature, the 
d x 2_ y 2 + d xy phase is created at temperature T = T c \. We study the temperature dependencies 
of the order parameter, specific heat and spin susceptibility in these mixed-angular-momentum 
states on square lattice and on a lattice with orthorhombic distortion. The above-mentioned phase 
transitions are identified by two jumps in specific heat at T c and T c \. 
PACS number(s): 74.20.Fg, 74.62.-c, 74.25.Bt 
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Inspite of many theoretical and experimental studies 
on high-T c cuprates the exact symmetry of the order pa- 
rameter is still a subject of active research [Q. However, 
there is evidence that the cuprates have singlet d-wave 
Cooper pairs and the order parameter has d x 2_ y i sym- 



metry in two dimensions |L| . Recent measurements M of 
penetration depth and superconducting specific heat at 
different temperatures T and related theoretical analyses 
HQ also support this. However, several phase-sensitive 
measurements of the order parameter of the cuprates 
indicate a significant mixing of a distinct angular mo- 
mentum component with a predominant d x 2_ y 2 state at 
temperatures below a second critical temperature T c \. 
For temperatures between T cl and T c only the d x 2_ y 2 
state survives. Below T c \ the order parameter can have a 
mixed-symmetry state of type d x 2_ y 2 + exp(i9)\, where 
X represents a state of different symmetry. The most 
probable possibilities for \ are the s or d xy wave. The 
possibility of a mixed (s — d)-wave symmetry was first 
suggested theoretically by Ruckenstein et al. and Kotliar 
1- 

There are experimental evidences based on Josephson 
supercurrent for tunneling between a conventional s-wave 
superconductor (Pb) and twinned or untwinned single 
crystals of YBa 2 Cu 3 07 (YBCO) that YBCO has mixed 
d x 2_ y 2 ± s or d x 2_ y 2 ± is symmetry || at lower temper- 
atures. Recently, the existence of these mixed-symmetry 
states has been explored to explain the nuclear magnetic 
resonance data in the superconductor YBCO and the 
Josephson critical current observed in YBCO-SNS and 
YBCO-Pb junctions Jjp. 

Kouznetsov et al. g performed some c-axis Josephson 
tunneling experiments by depositing conventional super- 
conductor (Pb) across a single twin boundary of a YBCO 
crystal. By measuring the critical current as a function 
of the angle and magnitude of a magnetic field applied in 
the plane of the junction they also found the evidence of 
a mixed-symmetry order parameter in YBCO involving 
d x 2_ y 2 and s waves. By measuring the microwave com- 
plex conductivity in the superconducting state of high 
quality YBa2Cu307_«5 single crystals at 10 GHz using 
a high-Q Nb cavity Sridhar etal also suggested the ex- 



istence of a multicomponent superconducting order pa- 
rameter in YBCO §. 

A similar conclusion of the existence of mixed- 
symmetry states may also be obtained based on the re- 
sults of angle-resolved photoemission spectroscopy exper- 
iment by Ma et al. in which a temperature dependent 
gap anisotropy in oxygen-annealed Bi2Sr2CaCu208+:r 
was detected [Hj]. The measured gaps along directions 
r — M and T — X are non-zero at low temperatures and 
their ratio was strongly temperature dependent. Us- 
ing Ginzburg-Landau theory, Betouras and Joynt Jll]] 
demonstrated that one way of explaining this behavior is 
to employ a mixed-symmetry state of the d x 2_ y 2 + s-wave 
type. They also conclude that the actual symmetry of the 
order parameter should vary substantially from one com- 
pound to another and for different levels of doping. This 
also suggests the possible appearance of a d x 2_ y 2 + d xy 
state under favorable conditions. 

More recently, Krishana et al. jl2) reported a phase 
transition in the superconductor Bi2Sr2CaCu20s in- 
duced by a magnetic field from a study of the thermal 
conductivity as a function of temperature and applied 
field. Laughlin |l3| provided a theoretical explanation 
of the observation by Krishana et al. || that for weak 
magnetic field a time-reversal symmetry breaking state 
of mixed symmetry is induced in Bi2Sr2CaCu20s. From 
a study of vortex in a d-wave superconductor using a 
self-consistent Bogoliubov-de Gennes formalism, Franz 
and Tesanovic |l3| also predicted the possibility of a su- 
perconducting state of mixed symmetry. This mixed- 
symmetry state is likely to be a minor s or d xy component 
superposed on a d x 2_ y 2 state for T < T c \. 

From different experimental observations it is now gen- 
erally accepted that a time-reversal symmetry breaking 
state of type d x 2_ y 2 + ix is possible in the presence of 
an external field or magnetic impurity. This mixed- 
symmetry state is observed close to these impurities, sur- 
faces/twin boundaries in the ab-plane or vortices. The 
nature of the mixed state varies from compound to com- 
pound. There are physical reasons for the appearance 
of these states. Either spin-orbit coupling with mag- 
netic impurities or Andreev reflected bound states which 
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create internal currents at the boundaries is responsible 
for these states JlJ]. However, orthorhombicity plays a 
crucial role in the generation of time-reversal symmet- 
ric mixed states. For example, it is established from 
a Ginzburg-Landau functional analysis Jig] that the or- 
thorhombicity has a consequence in the development of 
a d + s state instead of a time-reversal symmetry broken 
one. Moreover, from a theoretical point of view, time- 
reversal symmetric states of type d x 2_ y 2 +\ are expected 
to be allowed depending on the orthorhombic distortion. 

There have been some studies [|l6| on the phase tran- 
sition to a d x 2_ y 2 + exp(z0)x phase from a d x 2_ y 2 phase 
with 9 = 7r/2 and \ — d xy or a s state. From theoretical 
considerations we find that there are two possibilities for 
the phase 9: or ir/2. For 9 — 0, we find numerically 
that there is no stable d x 2_ y 2 + s phase. Here we study 
the phase transition to a d x 2_ y 2+d xy phase from a d x 2_ y 2 
phase below T c \ . In particular we study the temperature 
dependencies of the order parameter, specific heat, and 
spin susceptibility in the mixed-symmetry state. 

There is no suitable microscopic theory for high-T c su- 
perconductors and there is controversy about a proper 
description of the normal state and the pairing mecha- 
nism for such materials [jjj. In the absence of a micro- 
scopic theory, a phenomenological tight-binding model in 
two dimensions with the proper lattice symmetry will be 
used ]l7j . This model has been successful in describing 
many properties of high-T c materials. 

We study the temperature dependencies of specific 
heat and susceptibility of a d x 2_ y 2 + d xy -wave supercon- 
ductor with a weaker d xy wave both on square lattice and 
on a lattice with orthorhombic distortion. The order pa- 
rameter of a d x 2_ y 2 +d xy -wave superconductor has nodes 
on the Fermi surface and changes sign across it, and con- 
sequently, its superconducting observables also exhibit 
power-law dependencies on temperature. On the other 
hand, the order parameters for the mixed d x 2_ y 2 + is 
and d x 2_ y 2 + id xy -w&ve states do not have a node on the 
Fermi surface and the corresponding observables have a 
exponential dependencies on temperature. In the present 
study on d x 2_ y i + d xy -wave states the specific heat ex- 
hibits two jumps at T = T c \ and T = T c , which clearly 
exhibits the phase transition at T&. 

In the present two-dimensional tight binding model the 
effective interaction Vk q for transition from a momentum 
q to k is taken to be separable, and is expanded in terms 
of some general basis functions mu, labelled by the in- 
dex i, as Vkq = -J^iViVikVici [@- The functions rjiu 
are associated with a one dimensional irreducible repre- 
sentation of the point group of square lattice C^ v and 
are appropriate generalizations of the circular harmonics 
incorporating the proper lattice symmetry. The effec- 
tive interaction after including the two appropriate basis 
functions for singlet pairing is taken as 

Vkq = -V^lk^lq " V^k^q, (0.1) 

where rjiq = (cosq x — (3cosq y ) corresponds to d x 2_ y 2 



symmetry, r]2q = sin q x sin q y corresponds to d xy sym- 
metry, and where [3—1 corresponds to a square lat- 
tice, and ^ 1 represents orthorhombic distortion. In 
this case the quasiparticle dispersion relation is given 
by £k = — 2t[cosk x + (3 cos k y — jcosk x cosk y ), where t 
and [3t are the nearest-neighbour hopping integrals along 
the in-plane a and b axes, respectively, and "ft/2 is the 
second-nearest-neighbour hopping integral. The energy 
ek is measured with respect to the Fermi surface. 
At a finite T, one has the following BCS equation 
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with^q = [(eq-EV^ + IAql 2 ] 172 , where E F is the Fermi 
temperature and ks the Boltzmann's constant. The or- 
der parameter has the following anisotropic form: 



A q = Alfyiq + CA 2 772q, 



(0.3) 



where C is a complex nu mber of un it m odulas |C| = 
1. If we substitute Eqs. ([[!]) and (|o]| into the BCS 
equation (|0.2| ), one can separate the resultant equation 
in its real and imaginary parts. The resultant equations 
only have solution for real Ai and A 2 , when the complex 
parameter C is either purely real or purely imaginary. 
The solution for purely imaginary C, e.g., C = i have 
been extensively studied in relation to mixed d x 2_ y 2 + 
is and d x 2_„2 + id xy states [|l6|. Here we consider the 



d xy state. Usin g the 



solutio n f or C = 1, for d x 2_ y 

form (fl|) of A q with C = 1 and potential dOJJ), Eq. 
.2|) becomes the following coupled set of BCS equations 
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(0.5) 



where both the interactions V\ and V2 are assumed to 
be energy-independent constants for |e q — Ep\ < ksTu 
and zero for |e q — Ep\ > kBTjj, where ksTp, is a purely 
mathematical cutoff introduced to eliminate the ultravi- 
olet divergence in the BCS equation and should be com- 
pared with the physically motivated Debye cutoff in the 
case of the conventional superconductors 



We solved the coupled set of equations (0^) and (0.5) 
numerically and calculated the gaps Ai and A 2 at var- 
ious temperatures for T < T c . We have performed 
calculations (1) on a perfect square lattice and (2) in 
the presence of an orthorhombic distortion with cut off 
k B T D = 0.02586 eV (T D = 300 K) in both cases. The pa- 
rameters for these two cases are the following: (1) Square 
lattice - (a) t = 0.2586 eV, = 1, 7 = 0, V2 = 8.5t, and 
Vi = 0.734, T c = 71 K, T cl = 28 K; (b) t = 0.2586 eV, 
P = 1, 7 = 0, V 2 = 9.0t, and V x = 0.734, T c = 71 K, T cl 
= 55 K; (2) Orthorombic distortion - (a) t = 0.2586 eV, 
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(3 = 0.95, and 7 = 0, V 2 = 8.35*, and V 1 = 0.97*, T c = 70 
K, T ci = 30 K; (b) * = 0.2586 eV, (3 = 0.95, and 7 = 0, 
V 2 = 8.7*, and Vi = 0.97*, T c = 70 K, T cl = 50 K. For 
a very weak d x 2_ y 2-w&ve (d xy -wave) coupling the only 
possible solution corresponds to Ai = (A 2 = 0). We 
have studied t he so lution onl y when a coupling is allowed 
between Eqs. (D.4) and (X£). 
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Fig. 1. The order parameters Ai, A2 in Kelvin at 
different temperatures for models 1(a) (full line) and 1(b) 
(dashed line) on square lattice described in text. 

Figure 2 
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Fig. 2. The same as in Fig. 1 at different temper- 
atures for models 2(a) (full line) and 2(b) (dashed line) 
in the presence of orthorhombic distortion described in 
text. 

In Figs. 1 and 2 we plot the temperature dependencies 
of different A's for the following two sets of d x 2_ y 2 + d xy 
wave corresponding to models 1 and 2, respectively. In 
all cases, with the lowering of temperature passed T c , the 
parameter Ai increases up to T = T c \. As T is lowered 
further, A2 becomes nonzero at T — T c \ and begins to 
increase. As temperature is lowered, both Ax and A 2 
first increase and then attain a constant value at zero 
temperature. 

The different superconducting and normal specific 
heats are plotted in Figs. 3 and 4 for square lattice [mod- 



els 1(a) and 1(b)] and orthorhombic distortion [models 
2(a) and 2(b)], respectively. In both cases the specific 
heat exhibits two jumps — one at T c and another at T c \. 
From Figs. 1 and 2 we see that the temperature 

Figure 3 
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Fig. 3. Specific heat ratio C(T)/C n (T c ) versus T/T c 
for models 1(a) and 1(b) on square lattice: 1(a) (full line), 
1(b) (dashed line), d xy (dotted line), normal (dashed- 
dotted line). 

Figure 4 
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Fig. 4. Specific heat ratio C(T)/C n (T c ) versus T/T c 
for models 2(a) and 2(b) in the presence of orthorhombic 
distortion: 2(a) (full line), 2(b) (dashed line), d xy (dotted 
line), normal (dashed-dotted line). 

derivative of |A q | 2 has discontinuities at T c and T c \ due 
to the vanishing of Ai and A2, respectively, responsible 
for the two jumps in specific heat (see, definition in Ref. 
Q). For a pure d x 2_ y 2 wave we find that the specific 
heat exhibits a power-law dependence on temperature. 
However, the exponent of this dependence varies with 
temperature. For small T the exponent is approximately 
2.5, and for large T (T -> T c ) it is nearly 2. For the 
mixed d x 2_ y 2 + d xy -w&ve model, for T c > T > T c \ the 
specific heat exhibits d-wave power-law behavior. For d- 
wave models C s (T c )/C n (T c ) is a function of T c and (3. In 
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Figs. 2 and 3 this ratio for the d x 2_ y 2-wsve case, for T c = 
70 K, is approximately 3 (2.44) for = 1 (0.95). For the 
d^y-wave case, for T c = 70 K, this ratio is approximately 
1.81 (1.9) for fj = 1 (0.95). In a continuum calculation 
this ratio was 2 in the absence of a van Hove singularity 
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Fig. 5. Susceptibility ratio x( T )/x( T c) versus T/T c 
for square lattice: pure d x 2_ y 2 wave (solid line), pure d xy 
wave (dashed line), model 1(a) (dotted line), and model 
1(b) (dashed-double-dotted line). 



Figure 6 
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Fig. 6. Susceptibility ratio x( T )/x( T c) versus T/T c 
in the presence of orthorhombic distortion: pure d x 2_ y 2 
wave (solid line), pure d xy wave (dashed line), model 
1(a) (dotted line), and model 1(b) (dashed-double-dotted 
line). 

Next we exhibit the temperature dependence of spin 
susceptibility (defined in Ref. Q) in Figs. 5 and 6 where 
we also plot the results for pure d x 2_ y 2 and d xy waves 
for comparison. In Figs. 5 and 6, we show the results for 
models 1 and 2 on square lattice and with orthorhombic 
distortion, respectively. For pure d x 2_ y 2 and d xy waves 
we obtain power-law dependencies on temperature. The 
exponent for this power-law scaling was independent of 
critical temperature T c but varied from a square lattice to 



that with an orthorhombic distortion. For d x 2_ y 2 wave, 
the exponent for square lattice (orthorhombic distortion, 
(3 = 0.95) is 2.6 (2.4). For d xy wave, the exponent for 
square lattice (orthorhombic distortion, (3 — 0.95) is 1.1 
(1.6). For the mixed d x 2_ y 2 + d xy wave these exponents 
are nearly identical to the pure d x 2_ y 2 wave case. Hence, 
by studying the temperature dependency of spin suscep- 
tibility, it will be impossible to detect the phase transition 



d^y-wave 



at T = T c i from a d x 2_ y 2 wave to a d x 2_ y 2 + d xy 
at least within the present tight-binding model. 

In conclusion, we have studied the d x 2_ y 
superconductivity employing the two-dimensional tight 
binding BCS model on square lattice and also on a lat- 
tice with orthorhombic distortion and confirmed a second 
second-order phase transition at T — T c \ in the presence 
of a weaker d xy wave. This phase transition is marked by 
a jump in the specific heat at T = T c ±. We have kept the 
s- and d-w&ve couplings in such a domain that a coupled 
d x 2_ y 2 + (i^y-wave solution is allowed. The d x 2_ y 2 + d xy - 
wave state is similar to a da2_y2-wave-type state with 
nodes on the Fermi surface in the order parameter. Con- 
sequently, we find power-law temperature dependencies 
of specific heat and spin susceptibility in the d x 2_ y 2 +d xy 
wave. The exponents of these power laws for the mixed 
d x 2_ y 2 + d xy wave are very close to those for the pure 
d x 2_ y 2 wave. 
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